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Neuronal membrane properties can largely vary even within distinct morphological cell 
classes. The mechanisms and functional consequences of this diversity, however, are 
little explored. In the medial superior olive (MSO), a brainstem nucleus that performs 
binaural coincidence detection, membrane properties at rest are largely governed by 
the hyperpolarization-activated inward current <y which enables the temporally precise 
integration of excitatory and inhibitory inputs. Here, we report that lh density varies along 
the putative tonotopic axis of the MSO with lh being largest in ventral, high-frequency 
(HF) processing neurons. Also lh half-maximal activation voltage and time constant are 
differentially distributed such that ^ of the putative HF processing neurons activate faster 
and at more depolarized levels. Intracellular application of saturating concentrations of 
cyclic AMP removed the regional difference in hyperpolarization-activated cyclic nucleotide 
gated (HCN) channel activation, but not lh density. Experimental data in conjunction with a 
computational model suggest that increased ^ levels are helpful in counteracting temporal 
summation of phase-locked inhibitory inputs which is particularly prominent in HF neurons. 
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INTRODUCTION 

Neuronal encoding of information in the time domain is 
enhanced by specific adjustments of membrane properties to the 
dynamics and temporal characteristics of the inputs (O'Donnell 
and Nolan, 2011). This is especially important for neurons in 
the medial superior olive (MSO), a binaural nucleus in the audi- 
tory brainstem that analyses interaural time differences (ITDs) of 
different input frequencies with extremely high temporal preci- 
sion. This acuity primarily relies on the coincidence detection of 
precisely timed excitatory inputs from both ears onto MSO neu- 
rons (Grothe et al., 2010). In addition, two glycinergic inputs, 
originating from the ipsilateral medial and lateral nucleus of the 
trapezoid body, provide a prominent and phase-locked inhibition 
to MSO neurons, which fine-tunes the slope of the ITD func- 
tion to occur within the physiological range (Brand et al., 2002; 
Pecka et al, 2008; Leibold, 2010). Equally important for the high 
temporal precision with which these neurons integrate their exci- 
tatory and inhibitory inputs are the large voltage-gated channels 
that are open around the resting potential of the membrane. Such 
exquisitely fine-tuned temporal processing crucially depends on 
the composition and the properties of voltage-gated ion chan- 
nels. One of these voltage-gated currents is lh (or HCN-current), 
a cationic current, which is activated upon hyperpolarization 
(Wahl-Schott and Biel, 2009). Ih is especially large in MSO neu- 
rons and is regulated by intrinsic modulators such as cAMP 
and PIP2 (Khurana et al., 2012). In addition, these neurons also 
express a large low voltage-activated K + -channel (Klva) that also 
opens around the resting potential (Barnes-Davies et al., 2004; 
Mathews et al., 2010; Khurana et al., 2011). The sophisticated 



interplay between these channels reduces the input resistance and 
shortens the membrane time constant and thereby enhances the 
temporal acuity with which these neurons integrate their synaptic 
inputs (Barnes-Davies et al, 2004; Hassfurth et al, 2009; Mathews 
et al., 2010; Karcz et al., 2011; Khurana et al, 2011). 

Like most nuclei in the auditory brainstem, the MSO is tono- 
topically organized: Low-frequency (LF) sounds are represented 
dorsally and higher frequencies are processed ventrally (Guinan 
et al., 1972; Miiller, 1990). This spatial gradient of input fre- 
quencies enabled us to investigate the relationship between Ih 
properties, the integration of inhibitory inputs and its depen- 
dence on input frequency in the acute brain slice preparation 
using whole-cell patch-clamp recordings. We found that Ih is dif- 
ferentially distributed along the dorsoventral axis of the nucleus 
and that this spatial arrangement is paralleled by differential 
properties of synaptic integration. 

Moreover, we explored the putative functional consequences 
of this relationship theoretically using a computational single- 
compartment model featuring HCN and Klva channels that was 
fitted to electrophysiological recordings: this model suggests that 
integration of inhibitory inputs in a frequency-dependent man- 
ner helps to maintain the neuron's membrane potential close to 
firing threshold. 

MATERIALS AND METHODS 

All experiments were performed in accordance with the rules 
laid down by the EC Council Directive (86/89/ECC) and 
German animal welfare legislation and approved by the Regierung 
Oberbayern (AZ 55.2-1-54-2531-57-05, Bavaria, Germany). 
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All agents were purchased from Sigma-Aldrich (Germany) and 
Biotrend (Germany) unless otherwise indicated. 

SLICE PREPARATION 

Patch-clamp recordings were performed from MSO neurons of 
gerbils (Meriones unguiculatus) at the age of postnatal day 17/18 
(denoted P18) and 21/22/23 (denoted P22). The animals were 
decapitated under isoflurane anesthesia. The brain was removed 
in ice-cold oxygenated (95% 02/5% CO2) sucrose replace- 
ment solution containing (inmM): 2.5 KC1, 1.25 NaH2P04, 26 
NaHC0 3 , 0.25 CaCl 2 , 3 MgCl 2 , 12.5 glucose and 100 sucrose 
(pH7.4). Transverse brainstem slices (180 u,m) comprising the 
MSO were cut with a vibratome (VT1200S; Leica, Germany), 
incubated at 32°C for 15min in oxygenated artificial cere- 
brospinal fluid (ACSF) containing (inmM) 125 NaCl, 2.5 KC1, 
1.25 NaH 2 P0 4 , 26 NaHC0 3 , 2 CaCl 2 , 1 MgCl 2 and 25 glu- 
cose and then maintained at room temperature. For recordings, 
slices were transferred to a recording chamber, which was per- 
fused continuously with oxygenated ACSF at 32° C, and visualized 
with an upright microscope (Axioscope, Zeiss, Germany) using 
infrared-differential interference contrast optics. 

ELECTR0PHYSI0L0GY 

Current- and voltage-clamp recordings were made from visually 
identified MSO cells using a Multiclamp 700 A amplifier (Axon 
Instruments, USA) with standard electrode solution containing 
(inmM): 125K-gluconate, 5 KC1, 10 HEPES, 1 EGTA, 2 Na 2 ATP, 
2 MgATP, 0.3 Na 2 GTP and 10 Na-phosphocreatine; adjusted to 
pH 7.25 with KOH. All experiments were performed at near- 
physiological temperature (32°C). Patch pipettes were pulled 
from borosilicate glass capillaries (BioMedical Instruments, 
Germany) on a DMZ Universal Puller (Zeitz Instruments, 
Germany). When filled with electrode solution, patch pipettes 
had a resistance of 2-4 M£2. In some experiments Alexa-488 
(100 (iM) (Molecular Probes, Germany) was added to the elec- 
trode solution in order to verify the location of the neuron along 
the presumed tonotopic axis. 

During voltage-clamp recordings, whole-cell capacitance was 
compensated and used as measure for cell surface. The series 
resistance (<10M£2) was compensated to a residual of 2-2.5MA 
and not allowed to change more than 20%. To isolate Ih phar- 
macologically we applied the following drugs (inmM): 1 3,4 
diaminopyridine, 10 TEA-C1, 0.2 BaCl 2 , 0.001 TTX, 0.05 NiCl 2 , 
0.1 CdCl 2 , 0.01 DNQX, 0.025 DL-AP5 and 0.001 strychnine. NaCl 
was reduced to maintain iso-osmolarity. 

We cannot exclude that our voltage-clamp recordings are dis- 
torted due to space-clamp errors which result in incomplete con- 
trol of dendritic membrane potential. We minimize these errors 
by using an Ih isolation cocktail. Additionally, MSO neurons 
are anatomically compact cells with short dendrites (~150 |xm) 
(Rautenberg et al., 2009) so that space-clamp errors should be 
small. Moreover, it is likely that the somatic voltage-clamp under- 
estimates the HCN channel conductance. 

During current-clamp experiments, the bridge-balance was 
adjusted to compensate for artifacts arising from electrode resis- 
tance. In some experiments, Ih was blocked with the HCN 
channel-selective inhibitor ZD7288 (20 (jlM). 



Synaptic currents were evoked stimulating the slice with a glass 
electrode filled with 2 M NaCl. Stimulation electrodes were placed 
medial and lateral to the MSO. Inhibitory postsynaptic currents 
(IPSCs) were isolated by addition of 10 u,M DNQX and 25 |xM 
DL-AP5. IPSCs were evoked by brief pulses (100 u,s, intensities 
10-40 V) triggered by an analogue stimulus isolation unit (BST 
950, Dagan Corporation, USA). Patch electrodes were filled with 
(in mM) 99 CsMeS0 4 , 41 CsCl, 10 HEPES, 10 EGTA, 2 Na 2 ATP, 
2 MgATP, 0.3 Na 2 GTP, 5 TEA-C1 and 1 CaCl 2 , and 5 QX314 
to block postsynaptic Na + channels; adjusted to pH 7.25 with 
CsOH. 

In conductance-clamp experiments, simulated inhibitory con- 
ductances at 100 Hz were injected into MSO neurons with 
a SM-1 amplifier (Cambridge Conductance, UK). The simu- 
lated inhibitory conductance based upon recorded IPSCs (decay 
time: ~1.5ms, 10-90% rise time: ~0.9ms, reversal poten- 
tial: — 90 mV). The reversal potential was chosen according to 
data by Magnusson et al. (2005). 

DATA ACQUISITION AND ANALYSIS 

Both voltage and current signals were low-pass filtered at 10 kHz 
with a four-pole Bessel filter and sampled at a rate of 20-50 kHz. 
Stimulus generation and recordings were done with pCLAMP 
(Axon Instruments, USA). All electrophysiological data were 
analysed in IGOR Pro (Wavemetrics, USA) using Neuromatic 
and custom-written routines, or in Clampfit (Axon Instruments, 
USA). A junction potential of — 10.5 mV was corrected. 

Steady-state current responses were evaluated at the end of the 
voltage pulse. Ih density was obtained by normalizing the ampli- 
tude to the compensated whole-cell capacitance. The voltage 
dependence of Ih activation was measured from the tail cur- 
rent. Values were fitted with a Boltzmann function to obtain the 
half-maximal activation voltage Vn.5:/(V) = 1/(1 + exp[(Vn.5 — 
V/k)]), where V is the membrane voltage and k is the slope 
factor. The membrane time constants were evaluated by fit- 
ting a double-exponential function to the current traces: / (r) = 
Ai exp(-f/t fest ) + A 2 exp(-t/t s iow) where Tf ast and t s i ow are the 
fast and slow time constant of Ih activation. The effective time 
constant of Ih activation, T we ; g ht e d; was calculated according to: 

^-weighted 

(Ai * t fest + A 2 * t slow )/(Ai + A 2 ). V0.5 and 

^-weighted 

were estimated for each experiment and averaged. 

Input resistance was assessed from the peak hyperpolarization 
triggered by — lOOpA current injection according to Ohm's law 
R = U/I. The membrane time constant was estimated from a 
single-exponential fit to the voltage response to — 100 pA current 
injection. 

To determine decay times of evoked IPSCs the decay was fitted 
with a single-exponential function. The time course of inhibitory 
postsynaptic potentials (IPSPs) was analyzed by averaging 30 
traces, normalizing the resulting trace to the first IPSP amplitude, 
and then the 10-90% rise time, the 90-10% decay time and the 
half-width of the IPSPs were estimated. 

Results are expressed as mean ± standard error of the mean 
(SEM). Statistical significance was determined by a single-factor 
ANOVA test followed by a Scheffe's post-hoc test or by Student's 
unpaired f-test in Excel (Microsoft) with significance thresholds 
of P < 0.05 (*), P < 0.01 (**), and P < 0.001 (***). 
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RECONSTRUCTION OF PATCHED NEURONS 

Following recording, slices were fixed in 4% paraformaldehyde 
for 30 min. After extensive washing in phosphate-buffered saline 
(PBS) slices were exposed to blocking buffer (0.5% trition X- 
100/0.1% saponin/1% BSA in PBS) followed by incubation with 
the primary antibody (chicken anti-microtubule-associated pro- 
tein 2, MAP2, 1:1000, Neuromics) in blocking buffer. Slices were 
then rinsed in washing buffer (0.5% Trition X-100/0.1% saponin 
in PBS) and immunoreactivity was visualized by incubating the 
slices with the Cy3-conjugated secondary antibody raised in don- 
key (1:300; Dianova). Finally, slices were washed and mounted 
on slides with vectashield mounting reagent (Vector Laboratories, 
USA). 

MODELING 

A Hodgkin-Huxley-type single-compartment model was imple- 
mented separately for prototypic P22 dorsal and ventral cells. 
The temporal evolution of membrane potential V followed the 
differential equation 

dV 

Cm —j~ = — (4 + ^KLT + ^syn + ^1 eak ) 

with membrane capacitance C m and Ohmic currents 

I x (V)=g x afb n x (,V-E x ). 

The parameter g x describes the peak conductance, a x and b x are 
the gating variables for activation and inactivation, respectively, 
and E x denotes the reversal potential. The gating variables follow 
first order kinetics 

da Aqo — a db b^ — b 

— = and — = 

at x a at Xb 

with the steady-state activation a/boo and the voltage-dependent 
time constants x a /b. 

The low-threshold potassium channel (KLT) was modeled 
according to Mathews et al. (2010) with Ejc = — 90 mV. The 
kinetics of the hyperpolarization-activated cation current (Ih) was 
fitted to the data of voltage-clamp experiments from Figure 5, 
which resulted in the steady-state activation and the activation 
time constant (see Figure 7A) 

aoo (V) = (l + eO-nv + WAiy 1 and Tfl = 79 + 417e -(V+61.5) 2 /800 
for dorsal cells, and 

flco (V0 = (l + e 0.095( V + 75.5)^- 1 and Ta = 65 + 292fr (V + 62.5)V722_ 

for ventral cells, respectively (V inmV). Since HCN channels do 
not spontaneously inactivate, b was set to 1. As reversal potential 
we used = —35 mV. 

The model has been adapted to the different mean values 
of the membrane properties of the ventral (HF) and dorsal 
(LF) population by using the following channel peak con- 
ductances (innS/ii-m 2 ): g KLT dorsai = 0.0531, g HCN dorsai = 0.01025, 



gjQjpventral — gj^pdorsal * 5.4 and g HCN ventral = g HCN dorsal * 3.15. 

These settings yield a resting potential of around —60 mV for both 
model types and input resistances of _R; n = 23.94 MQ, for the dor- 
sal and R m = 3.77 M£2 for the ventral model corresponding to 
membrane time constants of x m of 1.64 ms and 0.45 ms, respec- 
tively. Using a specific membrane capacitance of 1 |xF/cm 2 these 
correspond to a modeled cell surface of 6839 |im 2 (dorsal) and 
12064 (xm 2 (ventral) with membrane capacities of 68.39 pF and 
120.64 pF, respectively. 

For both cell models the passive leak conductance was set to 
gleak = 33.3 fS/|xm 2 and the reversal potential was set to —70 mV. 

The fitting procedure described above implicates that the 
model parameters (specifically the HCN conductances) are 
adjusted according to our current-clamp data. This was done 
on purpose, since we assume the voltage-clamp data to be less 
accurate due to the above mentioned incomplete voltage-clamp 
control especially in the dendrites. 

The inhibitory input to the model was implemented as a 
conductance with reversal potential of —90 mV. The IPSG kinet- 
ics were fitted with a double-exponential (r in ms) to resemble 
measurements from Couchman et al. (2010): 

(l _ e -t/0.4) g-t/1.6 

G(t) = ginh max((l- e - f /°- 4 ) e - t / 1 - 6 ) 

For simulations to investigate the IPSP half-widths inhibitory 
100 Hz input stimuli of 20.5 nS (dorsal) and 90 nS (ventral) were 
applied to the models to roughly fit the membrane potential 
deflection seen in the corresponding current clamp experiments. 
The stimulus train was kept up for 800 ms to show the influence 
of the slowly activating HCN current. 

RESULTS 

l H VARIES ALONG THE D0RS0VENTRAL AXIS IN THE MSO 

Neuronal processing in the auditory system is tonotopically 
organized such that frequencies are orderly represented across 
most auditory nuclei. In the MSO, low frequency sounds are 
supposed to be encoded in the dorsal part of the MSO and 
higher frequency sounds are presumably represented in the 
ventral part (Guinan et al., 1972). In general, best frequen- 
cies of MSO neurons are lower compared to neurons in the 
LSO, but can occasionally still be above 2 kHz (Pecka et al., 
2008). Here, we investigated in a brain slice preparation of P18 
gerbils the biophysical properties of MSO neurons along this 
putative tonotopic axis. The MSO was subdivided into three 
regions, a ventral region, which we refer to as high-frequency 
(HF), a dorsal region, which we refer to as low-frequency (LF) 
and an intermediate middle-frequency (MF) region. MSO neu- 
rons were identified on the basis of their bipolar shape and 
their arrangement in a parasagittal plane. In some experiments, 
100 \iM Alexa-488 was included in the pipette solution to ver- 
ify the visually determined location of the neurons along the 
dorsoventral axis (Figure 1A). The properties of Ih between the 
regions were analyzed using voltage-clamp experiments. In all 
cells hyperpolarizing voltage pulses triggered slowly activating, 
large inward currents. Ih amplitude was 57% larger in ventral 
(presumably HF) neurons compared with dorsal (presumably 



Frontiers in Neural Circuits 



www.frontiersin.org 



July 2013 | Volume 7 | Article 117 | 3 



Baumann et al. 



Tonotopic organization of It, in the MSO 



dorsal J20MHL 



ventral 



dorsal 



ci rr 




< 

a 



c 
■a 



C2 



SO' 
_ 0- 
~ -50 • 
-100- 
-150- 
-200 

200 

2 £ 100 
t i 

3 Q. 50 



. #• o 

I** ° 
I** ° 



O ventral 

intermediate 
• dorsal 



-60 -40 mV 



D1 



300 
250- 
200 
150- 
100- 
50 
0 







m 


12 




18 









15 




12 



T 



E1 



, dorsal 




ventral intermediate dorsal 

ventral 



18 




D2 



ventral intermediate dorsal 



6OO-1 
500 
400 
300 
200 
100 
0 



E2 



I i 



i i 

n o " 



-130 -120 -110 -100 




E3 



-65 



ventral intermediate dorsal 



_ ^ -70 - 
If , 

s| ■' 

E « -f 



5 -85- 
-90- 



1 



FIGURE 1 | lh varies systematically along the dorsoventral axis. 

(A) A brain slice containing the MSO with Alexa-488-filled neurons (green) 
verifies the distribution of the patched neurons along the dorsoventral axis 
(red: MAP-2). (B) Pharmacologically isolated lh current traces were elicited 
by depolarizing and hyperpolarizing voltage steps from — 60.5mV to 
potentials between — 40.5mV and — 120. 5mV for 1 s in 5mV step 
increment and then to — 100. 5mV for 0.5s to elicit the tail current to 
determine the voltage dependence of lh activation. Current traces are 
representative for the dorsal, the intermediate and the ventral part of the 
MSO. (C) l-V relationships of steady-state (red arrow in B) l h density for 
ventral (n= 15), intermediate (n= 12) and dorsal (n= 18) neurons 
emphasize that lh density amplitudes are smallest in dorsal neurons and 



largest in ventral neurons (CI). Ij, density amplitudes for a voltage step 
to -110. 5mV (C2). (D) Weighted activation time constants at -110.5 mV 
(D1). The weighted activation time constants are voltage dependent and 
largest in the dorsal part of the MSO (D2). (E) The voltage-dependence of 
l h activation was measured from the tail current 20 ms after the end of the 
voltage steps (red arrow) (El). Values were fitted with a Boltzmann 
function to obtain the half-maximal activation voltage. In dorsal neurons 
the ^ activation curve is shifted to more negative voltages (E2). 
Half-maximal activation voltage was measured in each experiment and 
averaged (E3). Black symbols: dorsal neurons; gray symbols: intermediate 
neurons; white symbols: ventral neurons. **P < 0.01, ***P < 0.001, 
single-factor ANOVA test followed by a Scheffe's post-hoc test. 



LF) neurons [at -110.5mV: ventral: -3006 ±165 pA; inter- 
mediated: -2388 ± 123 pA; dorsal: -1910 ± 278 pA; ANOVA: 
F(2, 42) = 6.60, P = 0.003; Figure IB]. Also, lh density was signif- 
icantly larger in the ventral part of the MSO compared with the 
dorsal and intermediate part [at — 110.5 mV: ventral: — 127.7 ± 
5.3 pA/pF; intermediate: -85.1 ±5.4 pA/pF; dorsal: -73.4 ± 
9.3 pA/pF; ANOVA: F (2 , 42 ) = 14.86, P < 0.001; Figure 1C]. 
Dorsal neurons not only exhibited the smallest 1^ amplitude 
but also lh that activated slowest. The calculated weighted acti- 
vation time constant was two-fold larger in dorsal neurons 
compared with ventral neurons [at — 110.5 mV: ventral: 95.6 ± 
7.0ms; intermediate: 117.9 ± 9.2ms; dorsal: 191.3 ± 28.1ms; 



ANOVA: F(i, 42) = 6.62, P = 0.003; Figure 1D1]. For all three 
regions, the weighted activation time constants were voltage- 
dependent with T we i g hted = 72 ± 5 ms at — 120.5 mV increasing 
to T. we ighted = 215 ± 17 ms at —90.5 mV in the ventral part of the 
MSO (Student's paired f-test: P < 0.001) and with T we ; g hted = 
152 ± 22 ms at — 120.5 mV increasing to t we ighted = 367 ± 40 ms 
at —90.5 mV in the dorsal part (Student's paired f-test: P < 
0.001) (Figure 1D2). Analyzing the amplitude of the tail cur- 
rent revealed that lh voltage dependence was negatively shifted 
in dorsal neurons compared with ventral and intermediate neu- 
rons (Figure IE). Consequently, the half- maximal activation 
voltage was most negative in dorsal neurons [ventral: — 79 ± 
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1 mV; intermediate: -76 ± 2mV; dorsal: -87 ± 2mV; ANOVA: 
F(2, 42) = 13.51, P < 0.001; Figure IE]. On average, our mea- 
surements are in line with recently published data (Khurana et al., 
2012). 

Taken together, we observed a large difference in Ih proper- 
ties between the ventral and the dorsal part of the MSO. Dorsal 
neurons exhibited smaller Ih amplitude, slower activation kinetics 
and more negative half-maximal activation voltage as compared 
to ventral neurons. 

cAMP MODULATION OF \ h DIFFERS ALONG THE D0RS0VENTRAL AXIS 

HCN channel properties depend largely on the intracellular 
concentration of cAMP. The extent by which cAMP is able 
to regulate the gating of HCN channels is determined by the 
HCN subunits (Wahl-Schott and Biel, 2009). HCN1, which is 



less sensitive to cAMP, is the main subunit in MSO neurons 
(Koch et al, 2004; Khurana et al, 2012). Nevertheless, cAMP 
modulates the gating of HCN channels in the MSO proba- 
bly due to a co-assembly of HCN1 and HCN4 to heteromeric 
HCN channels (Khurana et al, 2012). To test whether a cAMP- 
dependent modulation underlies the differences in Ih proper- 
ties across the dorsoventral axis, we included 25 |xM cAMP in 
the pipette solution, which induces maximal cAMP modula- 
tion (Ludwig et al, 1998). As expected, Ih density amplitude 
was still significantly larger in ventral neurons compared with 
dorsal neurons (at -110.5 mV: ventral: —99.9 ± 6.6 pA/pF; dor- 
sal: -74.6 ± 6.8 pA/pF; Student's unpaired f-test, P = 0.014; 
Figures 2A,B). Moreover, cAMP accelerated the activation kinet- 
ics (Figures 2C,E) and positively shifted the activation curves in 
the two regions such that the activation curves overlapped for all 
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FIGURE 2 | Modulation of l h by cAMP differs along the dorsoventral axis. 

(A) Current responses to depolarizing and hyperpolarizing voltage steps were 
recorded from MSO neurons in the ventral and dorsal part of the MSO with 
25 |iM cAMP in the pipette solution. (B) l h density gradient persists in the 
presence of cAMP as illustrated by the current-voltage relationships for 
ventral (n = 13) and dorsal (n = 11 ) neurons (B1) and their l h density 
amplitudes for a —110.5 mV voltage step (B2). (C) The weighted activation 
time constants and (D), the voltage dependence of Ih activation overlap in the 



presence of cAMP Comparison of (E) the weighted activation time constants 
and (F) the half-maximal activation voltages in the absence and presence of 
25 (i,M cAMP reveals that dorsal neurons are more sensitive to cAMP than 
ventral neurons (F2). (F1) In the upper panel, tail currents were elicited using 
standard pipette solution. In the lower panel, a different, dorsal neuron is 
illustrated using standard pipette solution supplemented with 25|i,M cAMP 
Black symbols: dorsal neurons; white symbols: ventral neurons. *P < 0.05, 
»**p < 0.001 , two-tailed, unpaired f-test. 
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neurons (Figures 2D,F), with the largest shift observed in dor- 
sal neurons. Here, Weighted decreased more than two-fold from 
191.3 ± 28.1 ms(n= 18) to 92.1 ± 18.8 ms (n = 11) at 110.5 mV 
(Student's unpaired f-test: P = 0.012, Figure 2E) and half- 
maximal activation voltage increased by 15 mV from — 87 ± 2 mV 
(n = 18) to -72 ± 3 mV (n = 11) (Student's unpaired f-test: P < 
0.001). In ventral neurons, no shift of T we ; g hted was observed 
and half-maximal activation voltage was only shifted by about 
9mV (from -79 ± lmV (n = 15) to -70 ± lmV (« = 13); 
Student's unpaired f-test: P < 0.001; Figure 2F2). Hence, in the 
presence of saturating concentrations of cAMP the Ih activation 
kinetics and the dependence of Ih activation are similar whereas 
the dorsoventral difference of Ih amplitude persists. We, there- 
fore, assume that the spatial arrangement of Ih density originates 
from differences in HCN channel density, whereas distinct basal 
intracellular cAMP levels cause the dorsoventral organization 
of the half-maximal activation voltage and the activation time 
constants. 



I h DIFFERENCES AFFECT MEMBRANE PROPERTIES 

At rest a fraction of HCN channels is open in the dorsal part 
(~9%) as well as in the ventral part (~15%) of the MSO 
(Figure 1E2). This is in accordance with studies showing that 
Ih plays a critical role in determining the membrane proper- 
ties in auditory brainstem neurons (Golding et al., 1995; Adam 
et al, 2001; Koch and Grothe, 2003; Golding and Oertel, 2012). 
To test whether the observed differences in Ih result in diverse 
membrane properties we applied depolarizing as well as hyper- 
polarizing current injections and recorded the voltage responses 
from 59 neurons. As previously reported, depolarization of the 
cells elicited a single action potential at the onset of the cur- 
rent injection, whereas hyperpolarization induced a depolariz- 
ing voltage sag, which can be attributed to the activation of 
HCN channels (Figure 3A, Magnusson et al., 2005; Scott et al., 
2005). Despite the different open probability of HCN chan- 
nels at rest, the resting potential was nearly identical in all 
parts of the MSO (ventral: -58.8 ± 0.3 mV; dorsal: -59.1 ± 
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FIGURE 3 | Gradient of Ih affects membrane properties. (A) MSO neurons 
in the dorsal and in the ventral part of the MSO fire a single spike at the 
onset of depolarizing current injections and exhibit a voltage sag during 
hyperpolarizing current steps. (B) Resting membrane potential, (C) peak input 
resistance, and (D) membrane time constant are not significantly different 
between ventral (n = 1 8) and dorsal (n = 25) neurons (D2). The membrane 
time constant was fitted by a single-exponential function as shown in (D1 ). 
(E) Blocking If, with 20|lM ZD7288, a specific Ih blocker, hyperpolarizes the 



cell, increases the input resistance and the membrane time constant. Same 
ventral neuron as in (A) but after treatment with 20 \M ZD7288. The 
differences in (F) the input resistance (light-gray bars), the resting potential 
(dark-gray bars) and (G) the membrane time constant before and after 
application of 20 \lM ZD7288 varies between ventral and dorsal neurons. The 
effects are more pronounced in ventral neurons. Black symbols: dorsal 
neurons; white symbols: ventral neurons. *P < 0.05, two-tailed unpaired 
Mest. 



Frontiers in Neural Circuits 



www.frontiersin.org 



July 2013 | Volume 7 | Article 117 | 6 



Baumann et al. 



Tonotopic organization of I), in the MSO 



0.6 mV; Student's unpaired f-test: P = 0.692; Figure 3B) indicat- 
ing compensatory gradient of outward currents. The peak input 
resistance and the membrane time constant did not differ sig- 
nificantly between the frequency regions, however, both showed 
clear trends. Ventral neurons tended to exhibit the lowest input 
resistance (at -lOOpA: ventral: 10.7 ± 1.8MA; dorsal: 18.9 ± 
5.5 MS; Student's unpaired f-test: P = 0.240; Figure 3C). The 
membrane time constants were determined by fitting a single 
exponential function to the voltage traces (Figure 3D1). Ventral 
neurons tended to display the smallest membrane time constant 
(at -lOOpA: ventral: 0.69 ± 0.09ms; dorsal: 1.23 ± 0.27ms; 
Student's unpaired f-test: P = 0.108; Figure 3D2). To solidify the 
observed trends, we repeated the experiments under bath applica- 
tion of 20 u,M ZD7288, which selectively inhibits HCN channels. 
In all neurons, irrespective of their location along the dorsoven- 
tral axis, HCN channel blockade hyperpolarized the membrane 
potential and increased the input resistance and the membrane 
time constant (Figure 3E). This difference in input resistance and 
membrane time constant between control condition and HCN 
channel blockade varied significantly between dorsal and ven- 
tral neurons (Figures 3F,G). Thus, the fractional contribution 
of Ih is significantly different between dorsal and ventral neu- 
rons. The effects of 20 [iM ZD7288 were more pronounced in 
the ventral part of the MSO (Figures 3F,G) demonstrating that 
Ih contribution to the membrane properties is larger in ventral 
neurons, and confirming that the distinct membrane properties 
along the dorsoventral axis can be attributed to the observed 
differences in Ijj. 

INTEGRATION OF SIMULATED INHIBITORY INPUTS VARIES ALONG THE 
DORSOVENTRAL AXIS 

Assuming the MSO receives inputs that are phase-locked to the 
fine structure of a sound, the temporal summation of IPSP 
should vary between the regions, being most prominent in ventral 
neurons that presumably receive HF inputs and least in dor- 
sal neurons that presumably receive LF inputs. This summation 
would lead to a stronger hyperpolarization in ventral neurons 
and thereby reduce their excitability. In this case the observed 
dorsoventral difference of Ih, which activates upon hyperpolar- 
ization, would compensate for the putatively increased hyperpo- 
larization. To test our hypothesis, we simulated inhibitory inputs 
at 100 Hz and recorded the voltage responses from neurons in the 
dorsal and ventral part of the MSO. The simulated inhibitory con- 
ductance, which was injected into MSO neurons, was based upon 
recorded IPSCs (decay time: ~1.5 ms, 10-90% rise time: ~0.9 ms, 
amplitude: ~2 nA ). We also confirmed that decay times of IPSCs 
did not differ significantly between the regions (ventral: 2.2 ± 
0.1 ms, n = 9; dorsal: 2.7 ± 0.3 ms, n = 12; Student's unpaired f- 
test: P = 0.215; Figure 4B). These results are in line with data by 
Magnusson et al. (2005) showing that the IPSCs decay with time 
constants of around 1.5-3 ms in P18 gerbils (Magnusson et al., 
2005). 

As expected, the membrane potential response to the simu- 
lated IPSC trains varied as a function of the neuron's location 
along the dorsoventral axis (Figure 4A). The amplitude of the 
evoked IPSPs was significantly larger in neurons of the dor- 
sal part compared with neurons of the ventral part (ventral: 



18.3 ± 0.8 mV, n = 13; dorsal: 21.1 ± 0.9 mV, n = 10; Student's 
unpaired f-test: P = 0.035; Figure 4C1), also indicative for a 
larger input resistance in dorsal cells. Moreover, ventral neu- 
rons showed less summation than dorsal neurons (IPSP2/IPSP1: 
ventral: 0.980 ± 0.001, n = 13; dorsal: 0.995 ± 0.006, n = 10; 
Student's unpaired f-test: P= 0.011; Figure 4C2). To facilitate 
comparison of the time course, IPSPs were amplitude-normalized 
(Figure 4D1, inset) illustrating that the time course of IPSPs 
changed along the dorsoventral axis. The half-width of the IPSPs 
was largest in the dorsal part of the MSO and became smaller 
in the ventral part (ventral: 3.76 ± 0.25 ms, n = 13; dorsal: 
5.63 ± 0.60 ms, n = 10; Student's unpaired f-test: P = 0.005; 
Figure 4D). There was no difference in 10-90% rise time of 
the first IPSP between the frequency regions (ventral: 1.10 ± 
0.01ms; dorsal: 1.12 ± 0.03 ms; Student's unpaired f-test: P = 
0.450; Figure 4E), but the 90-10% decay time of the last IPSP was 
smallest in the ventral part (ventral: 2.51 ± 0.16 ms, n = 13; dor- 
sal: 4.05 ± 0.56 ms, n = 10; Student's unpaired f-test: P = 0.008; 
Figure 4E). Taken together, the time course of the IPSPs is faster 
in ventral neurons than in dorsal neurons. We speculate that these 
effects can be attributed to the dorsoventral organization of Ih 
as HCN channels are the main channel subtypes that open upon 
hyperpolarization. 

DORSOVENTRAL ORGANIZATION IS PRESERVED IN MORE MATURE 
ANIMALS 

There is evidence that HCN channels in the superior olivary 
complex undergo drastic developmental changes during the first 
three postnatal weeks (Leao et al., 2006; Hassfurth et al., 2009; 
Khurana et al., 2012). To rule out that these developmental 
refinements have implications on the observed dorsoventral orga- 
nization, we repeated experiments at P22, which is more at the 
end of this developmental period. Compared with P18, P22 MSO 
neurons exhibited slightly increased Ih amplitudes, the activa- 
tion curves were shifted to more positive half-maximal activation 
voltages and the activation kinetics were accelerated (Table 1). 
Nevertheless, Ih still varied systematically along the dorsoventral 
axis, such that ventral neurons exhibited significantly larger Ih 
amplitudes than dorsal neurons (at — 1 10.5 mV: ventral: —3304 ± 
1.5 pA, n = 8; dorsal: —2560 ± 238 pA, n = 9; Student's unpaired 
f-test: P = 0.019; Figure 5A). Accordingly, Ih density also var- 
ied significantly (at -110.5 mV: ventral; —122 A ± 9.6 pA/pF, n = 
8; dorsal: —92.7 ± 5.0 pA/pF, n = 9; Student's unpaired f-test: 
P = 0.012; Figure 5B). Also, the activation kinetics and the half- 
maximal activation voltage differed between the ventral part and 
the dorsal part of the MSO (tweighted at — 110.5 mV: ventral: 
78.0 ± 3.6 ms, n = 8; dorsal: 109.6 ± 18.7 ms, n = 9; Student's 
unpaired f-test: P = 0.138; Figure 5C; Vn.s: ventral: —76 ± 3 mV, 
n = 8; dorsal: — 81 ± 2 mV, n = 9; Student's unpaired f-test: P = 
0.160; Figure 5D). In P22 animals, Ih still was organized along the 
dorsoventral axis. However, these differences between ventral and 
dorsal neurons were less pronounced as compared to P18 gerbils. 

To assess to what extent these subtle changes in HCN proper- 
ties between P18 and P22 neurons affect the neurons' membrane 
properties, we also measured voltage changes in response to cur- 
rent injections in P22 animals. During depolarization neurons in 
the ventral part of the MSO fired a single spike at the beginning 
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FIGURE 4 | The integration of synaptic inputs varies along the 
dorsoventral axis. (A) Representative voltage traces to simulated IPSC 
trains (100 Hz). (B) The time course of IPSCs does not vary along the 
dorsoventral axis as indicated by the decay time. IPSCs were evoked by 
stimulating the slice medial or lateral the MSO with a stimulation 
electrode. (C) Dorsal neurons exhibit the largest IPSP amplitude (C1), 
and IPSP summation is increased in dorsal neurons (C2). (D) IPSP time 
course changes along the dorsoventral axis as depicted in the inset. The 



half-width of IPSPs is largest in dorsal neurons and the time course 
accelerates during stimulation with IPSC trains (D1). Half-width for the 
first IPSP (D2). (E) There is no difference between the rise times 
(10-90%) of the first IPSP (dark-gray bars) but the mean values for the 
decay time (90-10%) of the last IPSP (light-gray bars) are larger in dorsal 
neurons compared to ventral neurons. Black symbols: dorsal neurons; 
white symbols: ventral neurons. *P<0.05, **P<0.01, two-tailed paired 
or unpaired f-test, as appropriate. 



of the current injection (Figure 6A). In most ventral neurons, 
only for strong hyperpolarizing current injections the depolariz- 
ing voltage sag was obvious which is due to the extremely large 1^. 
This is also reflected in the very low input resistance (at — 100 pA: 
ventral: 3.7 ± 0.7 MS, n = 10; dorsal: 24.0 ± 6.4 MS, n = 12; 
Student's unpaired f-test: P = 0.016; Figure 6B) and in the very 
small time constant of ventral neurons (at — 100 pA: ventral: 
0.45 ± 0.07 ms, n = 10; dorsal: 1.64 ± 0.45, n = 12; Student's 
unpaired f-test: P = 0.047; Figure 6C). Compared with P18 ger- 
bils, the differences in the membrane time constant and in the 
input resistance between ventral and dorsal neurons were larger 
which resulted in significant differences along the dorsoventral 
axis (Figures 6B,C; Table 1). 

We evaluated the integration of inhibitory postsynaptic inputs 
by injecting currents with stimulus amplitudes adjusted to evoke 



physiological IPSPs of similar sizes (—8.1 ± 0.3 mV in the ven- 
tral part, n = 10, and —8.4 ± 0.3 mV in the dorsal part of the 
MSO, n = 12; Student's unpaired f-test: P = 0.493). Similar to 
PI 8, the voltage response to the simulated IPSC trains varied 
along the dorsoventral axis. The half-width of the first IPSP (ven- 
tral: 2.73 ± 0.06 ms, n= 10; dorsal: 4.29 ± 0.51ms, n = 13; 
Student's unpaired f-test: P = 0.014; Figure 6D), the 10-90% 
rise time of the first IPSP (ventral: 1.14 ± 0.03 ms, n = 10; dor- 
sal: 1.37 ± 0.07 ms, n = 6; Student's unpaired f-test: P = 0.007; 
Figure 6E) as well as the 90-10% decay time of the last IPSP 
(ventral: 1.84 ± 0.08 ms, n = 10; dorsal: 3.08 ± 0.39 ms, n = 13; 
P = 0.013; Figure 6E) were largest in the dorsal part of the MSO 
and became smaller in the ventral part. By comparing the time 
course of P18 and P22 neurons (Figure 6F, example for ventral 
neurons) we can demonstrate that consistent with an increase 
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Table 1 | Summary of HCN channel properties, membrane properties and synaptic properties of dorsal and ventral neurons for P18 and P22. 







P18 






P22 






Ventral 


Dorsal 


P 


Ventral 


Dorsal 


P 






Current [pA] 


-3006±165(15) 


-1909±277 (18) 




-3303 ±134 (8) 


-2559 ±238 (9) 


* 


Current density [pA/pF] 


-127.7±5.3(15) 


-73.4±9.3(18) 




122.4±4.9(6) 


92.7 ±5.0 (9) 


* 


Half-maximal activation voltage [mV] 


-79 ±1 (15) 


-87 ±2 (18) 




-76 ± 3 (8) 


-81 ±2 (9) 


n.s. 


T weighted I m s] 


95.6±7.0(15) 


191.3±28.1 (18) 


** 


78.0 ±3.6 (8) 


109.6 ±18.7 (9) 


n.s. 


MEMBRANE PROPERTIES (AT -100 pA) 


Resting potential [mV] 


-58.8 ±0.3 (24) 


-59.1 ±0.6(35) 


n.s. 


-60.0 ±0.7 (10) 


-60.1 ±0.9(12) 


n.s. 


Input resistance [MQ] 


10.7±1.8(24) 


18.9 ±5.5 (35) 


n.s. 


3.7 ±7.0 (10) 


24.0±6.4(12) 


* 


Membrane time constant [ms] 


0.69±0.09 (24) 


1.23 ±0.27 (35) 


n.s. 


0.45 ±0.07 (10) 


1.64 ±0.45 (12) 


* 



SYNAPTIC PROPERTIES 

IPSP half-width [ms] 3.76±0.25 (13) 5.63±0.60 (10) ** 2.72±0.06 (10) 4.29 ±0.51 (13) 

Rise time 10-90% [ms] 1.10±0.01 (13) 1.12±0.03 (10) n.s. 1 .14± 0.03 (10) 1.37±0.07 (13) 

Decay time 90-10% [ms] 2.51 ± 0.16 (13) 4.05± 0.56 (10) ** 1 .84± 0.08 (10) 3.08±0.39 (13) 

Values are mean ± SEM. The level of significance between ventral and dorsal was determined by using Student's unpaired t-test except for HCN channel properties 
of P18, where Scheffe's post-hoc test was employed following a single-factor ANOVA test (*P < 0.05, **P < 0.07, ***p < 0.007, n.s. not significant). The n-values 
for each group are stated in brackets. 



in Ih the half-width of the IPSPs is decreased in P22 neurons 
(Figure 6G). This emphasizes our hypothesis that Ih accelerates 
the time course of the IPSP and thereby decreases the temporal 
summation of IPSP. In addition, Ih compensates the summated 
hyperpolarization induced by the temporal summation of HF 
inhibitory inputs. 

Taken together, these data provide evidence that also in mature 
animals the integration of synaptic inputs varies as a function 
of the neuron's location along the dorsoventral axis and that a 
tonotopic organization of Ih may at least partially account for the 
observed gradient in synaptic integration. 

TONOTOPIC ORGANIZATION OF \ h ACCOUNTS FOR THE 
DORSOVENTRAL DIFFERENCES IN SYNAPTIC INTEGRATION 

To gain further mechanistic understanding and to assess the 
functional consequences of the dorsoventral Ih gradient in a 
computational model of a MSO cell, we first fitted activation 
profiles and channel time constants of Ih (from Figure 5) as 
described in the Materials and Methods section (Figure 7A). In 
addition to Ih, the model also included a low- voltage activated 
potassium current Ik-lva to counteract Ih induced depolariza- 
tion (Svirskis et al., 2002; see Materials and Methods). The peak 
conductances of Ih and Ik-lva were used as free parameters to 
adjust the neuron models to a given input resistance and rest- 
ing potential. Whereas the former was taken to be 3.77 M£2 for 
ventral (putative HF) neurons and 23.94 M£2 for dorsal (putative 
LF) neurons, the latter was assumed identical (— 60 mV) in both 
populations. 

We first validated our models by reproducing the current 
clamp experiments from Figure 6F (Figures 7B,C). Applying a 
stimulus of 100 Hz, we measured the half- width of the inhibitory 
potentials for cell models with both dorsal and ventral character- 
istics. The simulated IPSP half-widths are in very good agreement 
with the experimental data. 



Following the idea that the dorsoventral differences paral- 
lel the tonotopic axis, we simulated the response of the model 
neuron to periodic inhibitory inputs with different frequencies 
(Figure 7D). The kinetics of the individual IPSGs was modeled 
to fit those measured experimentally (Couchman et al., 2010 
and Materials and Methods). The decay constant x = 1.6 ms of 
these IPSGs is so slow that there is temporal summation of 
the IPSPs, which produces a significant hyperpolarizing volt- 
age offset (dark lines in Figure 7D). This offset increased for 
higher stimulation frequencies (200 Hz vs. 600 Hz in the exam- 
ple of Figure 7D). The increase of the hyperpolarizing volt- 
age offset for HF inputs can, however, be mitigated, if we 
assume that only the ventral neurons process HF inputs: In 
those neurons this offset is smaller because of the lower input 
resistance that results from larger Ih and Ik-lva conductances 
(Figure 7D). 

The kinetics of Ih are much slower than the time constants 
that are typical for fast auditory processing. Therefore, Ih is gen- 
erally assumed not to be suited to directly interact with neuronal 
processing of sound information on a fast time scale. However, 
the interplay between Ih and Ik-lva may play an important role 
in temporal sharpening of the PSPs (Khurana et al., 2011). In 
contrast to Ih, the Ik-lva does possess fast kinetics and thus 
has been proposed to contribute to fast temporal processing of 
MSO neurons in several studies (Svirskis et al., 2002; Jercog et al., 
2010). To specifically evaluate the interaction between Ih and 
Ik-lva channel kinetics in the present context of IPSG trains, 
we also simulated a model in which the Ik-lva kinetics had 
been slowed down such that the kinetics were comparable to 
the Ih kinetics, while leaving the input resistance unchanged. 
Comparing both models (fast Ik-lva kinetics, Figure 7C and slow 
Ik-lva kinetics, Figure 7E) we found a clear effect on temporal 
precision as measured by an increase in IPSP half-width. This 
increase was stronger for the model of the ventral MSO neuron 
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FIGURE 5 1 lh gradient persists in more mature animals (P22). 

(A) Pharmacologically isolated lh current traces were elicited by 
depolarizing and hyperpolarizing voltage steps from — 60.5mV to 
potentials between — 40.5mV and — 120. 5mV (5mV step increment). 
Current traces are representative for the dorsal and the ventral part 
of the MSO. (B) l-V relationships of steady-state If, density for 
ventral (n=8) and dorsal (n=10) neurons emphasize that lh density 
amplitudes are smallest in dorsal neurons and largest in ventral 



neurons (B1). lh density amplitudes for a voltage step to — 110. 5mV 
(B2). (C) The weighted activation time constants are voltage 
dependent and largest in the dorsal part of the MSO. (D) The 
voltage-dependence of lh activation was measured from the tail 
current. In dorsal neurons the lh activation curve is shifted to more 
negative voltages (D1). Half-maximal activation voltage was measured 
in each experiment and averaged (D2). Black symbols: dorsal neurons; 
white symbols: ventral neurons. *P<0.05, two-tailed unpaired f-test. 



(Figure 7E), i.e., a putative HF processing neuron, although the 
input resistance was the same for both Ik-lva kinetics. This shows 
that for high frequencies neurons with large lh the temporal pre- 
cision of the hyperpolarizing IPSPs is considerably enhanced by 
the active properties of the fast Klva channels, whereas for low 
frequency neurons this temporal integration is mostly explained 
by the differences in input resistance. Mechanistically, the In- 
dependent sharpening of IPSPs can be understood as follows: 
During the hyperpolarizing flank of the IPSPs the Ik-lva chan- 
nels - which are open at rest - close very rapidly and thereby 
effectively set a new equilibrium potential of the whole cell at 
a depolarized level close to the reversal of lh- The resulting 
huge driving force massively speeds up the depolarizing flank of 
the IPSP and thereby accounts for the temporal sharpening. As 
the membrane potential approaches the old equilibrium poten- 
tial, the Klva channels quickly open again and they restore the 
original equilibrium potential with only little overshoot as wit- 
nessed by the small amplitude of the voltage fluctuation after 
the IPSPs in Figure 7C. To test the above hypothesis, we con- 
ducted simulations with different reversal potentials of I^, As 



expected, a reduction of the driving force broadened the IPSPs 
(Figures 7F,G). 

In summary, we conclude that fast Klva channels in interac- 
tion with lh may predominantly sharpen the IPSPs (particularly 
in HF neurons with large 1^), whereas 1^ in MSO neurons 
alone balances out the hyperpolarizing voltage offset induced 
by the temporal summation of phase-locked inhibitory synaptic 
currents. 

DISCUSSION 

In the present study we demonstrate that 1^ amplitude systemati- 
cally varies along the dorsoventral axis of the MSO, being largest 
in ventral neurons and smallest in dorsal neurons. Consistent 
with this dorsoventral organization of membrane properties, the 
integration of inhibitory inputs systematically varies as a func- 
tion of the neuron's location in both experiments and the model 
indicating that MSO neurons are tuned differentially along the 
presumed tonotopic axis. Tonotopic gradients of lh have been 
previously observed in auditory brainstem nuclei. For example, 
in the lateral superior olive (LSO) lh is larger in the LF region of 
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FIGURE 6 | Differences in membrane properties are more evident in 
mature animals. (A) MSO neurons in the dorsal and in the ventral part of the 
MSO fire a single spike at the onset of depolarizing current injections and 
exhibit a voltage sag during hyperpolarizing current steps. In particular, in 
ventral neurons the voltage sag only is obvious for strong hyperpolarizing 
current injections due to the large resting conductance of If,. (B) Peak input 
resistance and (C) membrane time constant vary significantly between 
ventral (n= 10) and dorsal (n= 12) neurons. Differences in the input 
resistance and in the membrane time constant between ventral and dorsal 
neurons are significantly larger in P22 animals compared with P18 animals. 



(D) IPSP time course changes along the dorsoventral axis as depicted in the 
inset. The half-width of IPSPs is largest in dorsal neurons and the time course 
accelerates during stimulation with IPSC trains. (E) The mean values for the 
rise time (10-90%) of the first IPSP (dark-gray bars) and the decay time 
(90-10%) of the last IPSP (light-gray bars) are larger in dorsal neurons 
compared to ventral neurons. (F) IPSP time course is accelerated in P22 
compared to P18. Representative normalized IPSP trains for the ventral part of 
the MSO. (G) Half-widths for the first IPSPs are decreased for P22 in both the 
ventral and the dorsal part of the MSO. Black symbols: dorsal neurons; white 
symbols: ventral neurons. *P < 0.05, **P < 0.01, two-tailed unpaired f-test. 



the nucleus compared to the HF region (Hassfurth et al., 2009). 
This opposite gradient might be due to the fact that in general 
LSO processes much higher input frequencies compared to MSO 
neurons (Sanes et al., 1989; Tolnai et al., 2008). The Ih gradient 
is also opposite in the nucleus laminaris (NL) (Yamada et al., 
2005), the bird's MSO analogue ITD processing stage. Whether 
this difference is due to the diverse function of inhibitory inputs 
in mammals and birds or to different ITD processing strategies in 
the two animal classes is not clear. However, these results suggest 
that tuning of biophysical membrane properties through differ- 
ential expression of HCN channels along the tonotopic axis in 
general optimizes the processing of different inputs frequencies 
(Kuba et al, 2005; Slee et al, 2010). 

In mammals, Ih (or HCN) channels can derive from 
four different genes (HCN 1-4) and assemble into homo- or 



heterotetramers with distinct electrophysiological properties in 
terms of their activation kinetics, their activation dependence, 
and their sensitivity to cAMP. In contrast, single-channel conduc- 
tance is very similar for the different HCN isoforms (Brandt et al., 
2009) and maximal Ih amplitude depends only very little on intra- 
cellular modulators (Ludwig et al., 1998; Wahl-Schott and Biel, 
2009). This suggests that the observed dorsoventral gradient of Ih 
density in the MSO most likely relies on differences in the num- 
ber of HCN channels and is independent of subunit variation. 
Our experiments also show that Ih activation kinetics accelerates 
and half-maximal activation voltage increases from the dorsal to 
the ventral part of the MSO. A distinct distribution of the iso- 
forms along the dorsoventral axis might provide an explanation 
for the differences in biophysical properties of Ih (Yamada et al., 
2005). MSO neurons mainly express HCN1 and HCN4 subunits 
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FIGURE 7 | Integration of IPSG trains, computational model. (A) 

Activation curve and time constant of lh from Figure 5 (symbols) were fitted 
for the ventral and dorsal population (solid lines). (B) A 100 Hz IPSG train 
(green) was applied to the cell models representative for the dorsal (blue) and 
ventral (red) population (see Materials and Methods). The individual inhibitory 
conductance was set to 20.5 nS and 90 nS for the dorsal and for the ventral 
cell models to obtain similar voltage amplitudes. (C) Half-width of the 
resulting IPSPs from (B). The first nine IPSPs were magnified on the right 
and overlaid with the respective K LVA steady-state activation (gray levels). 



(D) Voltage responses (light colors) for different stimulus frequencies 
(as indicated) and model cells (colors). The dark traces were generated by 
low-pass filtering the voltage trace with a second-order Butterworth low-pass 
filter with cut-off frequency of 100 Hz. (E) Same as (C) for a Klva model with 
100-fold slowed down activation and inactivation time constants. 
Conductances for If, and Ik-lva were adjusted to match the input resistances 
and equilibrium potential of the cells in (C). (F) IPSP trains for different 
reversal potentials of ^ (color code see G). (G) Half-widths of the IPSP trains 
from (F). 



which both possess distinct physiological properties (Khurana 
et al, 2012). Among all different HCN subunits, HCN4 pos- 
sesses the slowest kinetics and the most negative half-maximal 
activation voltage, whereas HCN1 possesses the fastest kinetics 
and the most positive half-maximal activation voltage (Santoro 



et al., 2000; Moosmang et al., 2001). Thus, an increased contri- 
bution of HCN1 towards ventral neurons could result in faster 
activation kinetics and more positive half-maximal activation 
voltage. Conversely, our data suggest that different basal levels of 
intracellular cAMP cause the observed dorsoventral gradient in Ih 
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properties. The gating of HCN channels in MSO neurons is very 
sensitive to cAMP, since most likely HCN1 and HCN4 isoforms 
co-assemble to form fast-activating but cAMP-sensitive HCN het- 
eromers (Khurana et al., 2012). We show that dialyzing neurons 
with a saturating cAMP concentration resulted in nearly identi- 
cal activation kinetics and half-maximal activation voltages in all 
MSO neurons. This opens the possibility that differential activa- 
tion or expression of receptors that modulate intracellular cAMP 
levels could modify Ih properties along the presumed tonotopic 
axis of the MSO (Yamada et al., 2005) and regulate processing of 
various input frequencies in an activity dependent manner. 

Functionally, Ih strongly influences basic membrane proper- 
ties such as resting potential, input resistance and membrane 
time constant of neurons. These properties determine cellular 
excitability and synaptic integration. More specifically, Ih depo- 
larizes the resting potential toward spike threshold, decreases 
the membrane time constant and lowers the input resistance at 
and below the resting potential, when the membrane potential is 
hyperpolarized in response to inhibitory inputs. Consistent with 
this idea we found that ventral neurons had a lower input resis- 
tance and a faster membrane time constant than dorsal neurons. 
In addition, postsynaptic integration of inhibitory inputs differed 
dependent on Ih amplitude and properties. 

But what are the functional implications for the diversity 
of Ih of different neuron types for information processing in 
a small network? This and the relation to Ih has been exten- 
sively studied in both the entorhinal cortex and the hippocampus 
where Ih properties and HCN channels are as well distributed 
along a dorsoventral gradient (Garden et al., 2008; Giocomo and 
Hasselmo, 2008; Marcelin et al, 2012a,b). In these structures, Ih 
has been hypothesized to contribute to the observed gradient in 
grid field spacing in the entorhinal cortex (Giocomo et al., 2011; 
Hussaini et al., 2011). This mostly relates to the fact that Ih accel- 
erates resonance frequency in those neurons. In these neurons Ih 
also tunes the membrane properties to the slow oscillatory activity 
of the inputs they receive, which is crucial for the specific function 
of these neurons. In auditory brainstem neurons and especially 
in the MSO, input frequencies (up to 1.5 kHz) are a magnitude 
higher than the activation and deactivation kinetics of Ih and thus 
an active contribution to temporal processing is unlikely. 

One possible explanation why Ih distribution is tonotopically 
organized is suggested by our model and our experimental data. 
Neurons in the MSO not only receive two precisely timed exci- 
tatory but also two prominent inhibitory inputs from the medial 



and lateral nucleus of the trapezoid body (Grothe et al., 2010), 
which are phase-locked to the fine-structure of the sound. Due 
to the relatively slow time constants of the inhibitory inputs 
(Magnusson et al., 2005; Couchman et al., 2010), the inhibi- 
tion summates and strongly hyperpolarizes the neuron. Since Ih 
is rapidly activated during hyperpolarization we propose that Ih 
reduces the integration of synaptic inputs during periods of pro- 
longed hyperpolarization. Indeed, both our experimental data 
and our model show that Ih decreases the temporal summa- 
tion of IPSPs by gradually activating and thereby opposing the 
summated hyperpolarization induced by the temporal summa- 
tion of HF inhibitory inputs. Functionally such a hyperpolarizing 
offset is problematic, since it effectively increases the spike thresh- 
old and thereby strongly reduces or even completely prohibits 
neuronal spiking in response to these input frequencies. The 
additional Ih activated in the ventral MSO region prevents this 
excessive hyperpolarization and keeps the neurons in an operat- 
ing regime for binaural coincidence detection. In MSO neurons 
both Ih and Ik-lva are open at rest (Khurana et al., 2011) and 
both contribute to the extremely low membrane time constants. 
The balance between the hyperpolarizing Ik-lva and the depo- 
larizing Ih determines the resting potential and together lowers 
the membrane time constants in both the hyperpolarizing and 
the depolarizing range. This decrease in time constant also in 
the depolarizing range would then improve coincidence detection 
of inputs thereby optimizing ITD analysis in these neurons. In 
addition, a higher expression level of Ih also indirectly enhances 
Ik-lva -induced sharpening of inhibitory synaptic potentials by 
modulating the speed of depolarization via the driving force of 
Ih (Figures 7C,E). Conversely, two recent studies suggest that 
increasing Ih in the MSO and the NL, the bird's analogue structure 
of the MSO, sharpens the time window for coincidence detec- 
tion also of excitatory inputs (Yamada et al., 2005; Khurana et al, 
2012). This all implies that MSO neurons that respond best to 
higher frequency sounds and have thus larger Ih should have 
sharper time windows for ITD detection compared to neurons 
responding best to low frequency sounds. This phenomenon can 
indeed be observed for ITD functions of MSO neurons that are 
tuned to different best frequencies (Yin and Chan, 1990; Brand 
et al., 2002). 
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